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1. Introduction 

Recently, there has been great interest in studying the properties of confined quantum 
systems. There are many reasons for this. The recent developments in nanotechnology 
have generated intensive research activity in modeling spatially confined quantum 
systems [1]-|2]. When an atom or a molecule is trapped inside any kind of microscopic 
cavity, or is placed in a high pressure environment, it experiences spatial confinement 
that affects its physical and chemical properties [3]-[l]. The concept of a confined 
quantum system goes back to the early work of Michels et al [1] who studied the 
properties of an atomic system under very high pressures. They suggested to replace 
the interaction of the atoms with surrounding atoms by a uniform pressure on a sphere 
within which the atom is considered to be en closed. This led them to consider the 
problem of hydrogen with modified external boundary conditions [2]. Since then, the 
confined hydrogen atom attracted widespread attention [2] -[33]. The approach to this 
confined problem is the same as that used to study a hydrogenic donor located at the 
center of a spherical GaAs — {Ga, Al)As quantum dot, a semiconductor device that 
confines electrons [13]-[l7j. Considerable interest in calculating the properties of donor 
states in a quantum dot has been renewed recently and a number of calculations on 
the bound states of a hydrogenic donor in a quantum dot have been reported (pTj. 
and the references therein). Many researchers have carried out accurate calculations 
of eigenvalues of the confined hydrogen atom using various techniques. Some of these 
are variational methods [18]- [27], finite element methods [28], and algebraic methods 
[29] . In the present work, we have calculated the energy eigenstates of the confined 
hydrogen atom for different quantum levels using the Asymptotic Iteration Method 
(AIM) [31]. In addition to calculating the energy eigenvalues, we also used AIM to 
compute the critical cage radius at which a state is no longer bound. Sommerfeld and 
Welker [5] gave a detailed investigation on the variation of the binding energy of Is state 
of hydrogen atom with respect to the sphere radius. They showed that there is critical 
value of the sphere radius Vc at which the binding energy is zero. This radius is now 
known as the 'critical cage radius' [30j. It has been found that, for r < r^ the energy of 
the system is positive, i.e. the electron exerts pressure on the walls of the sphere [30j . 
Later some workers calculated this critical cage radius for various atoms and different 
quantum levels. For example, Boeyens [32] has calculated the Tc values for the ground 
state of many atoms by the Hartree-Fock-Slater method. 

The organization of the paper is as follows. In the next section, we give a brief 
outline of the confined Hydrogen-like atom problem. In section 3, we summarize the 
asymptotic iteration method (AIM) used in this work. In section 4, we apply AIM to 
the confined Hydrogen-like atom problem and derive certain exact solutions for special 
parameter values. In section 5, we present the results of our general calculations for 
any box radius R, and we compare them with the most accurate results available in the 
literature. We present our calculated critical cage radii in section 6. 



ip{r) = Ei){r), V^(0) = i){R) = (1) 
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2. Formulation of the Problem 

In atomic units, the radial Schrodinger equation of a Hydrogen-like atom located in the 
center of a spherical box of radius R can be written as 

_(P_ l{l + l) A' 

dr"^ ^2 J. 

where A > is real parameter, E is the atom's energy and the boundary condition 
for the exact eigenfunction \E'(r) is ^E'(O) = "^{R) = (i.e. it satisfies the Dirichlet 
boundary condition). Here, r is the electron position vector with respect to the nucleus 
and r = ||r|| is the length of this vector. Of course, the presence of the Dirichlet 
boundary condition changes the structure of the atomic energy spectrum drastically as 
compared with that of the free, non-confined, hydrogen atom corresponding to the limit 
R ^ oo. The interesting hidden symmetry that manifests itself in extra degeneracy of 
energy levels at specific values of R has been studied [S]-[I2] in number of articles. 
It is interesting to note that when r tends to zero, the wavefunction \l/(r) behaves as 
r'"*"^, and when r approaches R, \l/ must approach zero. Under these assumptions, we 
may represent the unnormalized wavefunction in the form 

, . J r'+^(_R — r) exp(— ar)/(r), ii r < R , , 

"^^""^^JO iir>R ^^^ 

where /(r) is to be determined, and a is a parameter that will be used to obtain the 
energy and to control the convergence of the iterative method used by AIM. As we shall 
show, this simple form of wave function ([21) gives good results for the confined hydrogen 
atom. Substituting ([2]) in the Schrodinger equation ([1]), we find that the radial function 
/(r) must satisfy the differential equation 

1 / + r 



f"{r) = 2 a + 



R — r r 




where we denote 

E = -al (4) 

The problem is then reduced to solving this second-order homogeneous differential 
equation for /(r). The asymptotic iteration method (AIM) was developed with the idea 
of using minimal algebraic computations to solve such differential equations. In the 
next section, we give a brief introduction to the asymptotic iteration method. Detailed 
proofs and applications to unconfined systems can be found in [31] and [35] . 

3. Brief Introduction to the Asymptotic Iteration Method 

Given Ao(a;) and so{x) sufficiently differentiable functions, the asymptotic iteration 
method tells us that the second-order differential equation 

y" = \o{x)y' + so{x) (5) 
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has a general solution 

y{x) = exp(- r adt) fe + Ci T exp ( f (Ao(r) + 2a(r))rfr') dt) (6) 



if for some n > 0, 

^n °n—l 



An A 



(7) 



n-l 



where 



^n = ■Sn-1 + SQAn-l- 



(8) 




In general, the (asymptotic) termination condition Eq.(7) can be written equivalently 
as follows 

Sn{r) = XnSn^l- SnXn-l = 0, 72 = 1,2,... (9) 

Note that, we can start the computation of the recurrence relation (Q from n = with 
the initial conditions A_i = 1 and s_i = 0. It follows that if Sn{r) = 0, then (5„+i(r) = 
for all n. The termination condition ([9]) has a crucial role for the computation of the 
eigenenergies (jl]). Indeed, using ([3]), we obtain, by means of 

= 2(a+ — ^±i) 

_ i (2l+2)a~A I 21+2 _ _2a_\ ^UJ 

V r ~'~ r(R-r) R-r > ^ 

the recursive relations for s„ and A„, n = 1, 2, . . . given by (jH])- These quantities are, in 
general, functions of the parameter a and the variable r. If, for a suitable choice of a, 
the termination condition Eq.Q is satisfied at every r < R, then the problem is called 
'exactly solvable'. If the differential equation is not exactly solvable, the function S{r), 
in general, will depends on both a and r. In this case, we obtain a by iterating Eq.Q 
with a suitable initial value of r = tq < -R. 

4. Exact analytical solutions 

In this section, we first present the exact analytical solutions for a Hydrogen-like atom 
confined in a spherical box of radius R. The application of AIM which we employ here 
shows that for certain values of the parameter a, we obtain exact solution provided the 
confinement radius R assumes corresponding definite values; that is to say, Eq.([3]) is 
then exactly solvable. We first note that the integer n here is not to be identified for 
the hydrogenic problem with the radial quantum number, but rather by n = l+number 
of radial nodes. Thus, with this notation, in the limit as R ^ oo, the eigenvalues of 
Eq.(l) become exactly E = —A^/{A{n + £)^). Using ([9]) and (ITOl) we find the following 
results. 

• For n = 1 and a = 20^5 'we have (5i = if 



AR - 2{l + l){l + 2) = 0. 
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A 
2(1+3) ■' 



For n = 2 and a = „,,^^ , we have ^2 = if 



A^R^ - 2(2/ + 3)(/ + 3)AR + 2{l + 3)^(2/ + 3)(/ + 1) = 0. 
For n = 3 and a = 2(i^a) i ^^ have ^3 = if 

A^R^ - 6(/ + 2)(/ + A)A^R^ + 6(/ + Af{l + 2)(2/ + ?,)AR 
-4(/ + 4)3(/ + 2)(/ + l)(2/ + 3) = 0. 



and so on 


. In general, we have for 


= 0,1,2,. 








"" 2{n + l + iy 




that 5n = 


0, n = l,2,... if 










^(-l)^(/ + n + lf- 


-*^r(2/ + 


n + 2)V{n 


+ 1) 



(11) 



^^ T{2l + k + 2)T{n + l-k)V{k + l) "^^^^ °" ^^^^ 

It is interesting to note that the polynomial conditions on R given by flT2l) can be written 
in terms of the confluent hypergeometric functions as 

iFi(-n;2/ + 2; /^ )=0, n = l,2,.... (13) 

V / + n + 1/ 

This means that, in order to obtain the energy spectrum as given by (jl]), we must find 

the roots of the confluent hypergeometric function (1131) . In this case, for given A, the 

eigenvalues (jl]) are given by 

A^ 
^"' = - ,, , , , n2 ' ^ = 1, 2, . . . , / = 0, 1, 2, . . . (14) 

4(n + f + 1)^ 

where R has values which are the roots of flT3|) . The corresponding analytic wave 
functions can be computed by using Eq.(l6]) for A and i?, related by means of (IT3|) for 
given n and /. Straightforward computations then show that 

fii{r) = l, fni{r) = Ti(l--].i^m (15) 

where m = 0,1,2, ...n — 1 and Ri are the roots of (fT3|) . In Table 1, we report some 
eigenenergies of the enclosed hydrogen atom as a function of the radius R along with the 
corresponding wave functions /(r) computed by using (IT5|l . We can obtain the complete 

wavefunction from Eq.([2]). For given n,l, m quantum numbers, the corresponding (un- 
normalized) wave functions are given by 



*- - ^- ^'" ->^ (-ndr^') n (^ - ^) • ^ ^ '»• «"' 



(16) 



where Ri{m) are given as the roots of the confluent hypergeometric function Eq. flT3l) . We 
have first obtained these exact solutions of the confined hydrogen-like problem, valid 
for special values of the parameters, since they can be very useful for verifying the 
correctness of general approximations. In the next section, we will give more general 
results for the problem by applying AIM to Eq. ^ for arbitrary given R. 
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Table 1. Exact energy eigenvalues for certain fixed values of R. 



n 


/ 


m 


R 


E 


/n7W 


1 








4 
A 


16 




1 


1 





12 
A 


A^ 
36 




1 


2 





24 
A 


64 




1 


3 





40 
A 


A^ 
100 




2 







1 


3(3- v^) 

A 
3(3+^3) 

A 


A' 
36 
A^ 
36 


1 Ar 




3(3+v/3) 
1 Ar 


3(3-v^) 


2 


1 




1 


4(5-^5) 

A 
4(5+^5) 

A 


^2 

64 
64 


1 Ar 




4(5+v^) 
1 Ar 


4(5-V5) 


2 


2 




1 


5(7- v^) 
A 

A 


A2 
100 

100 


1 Ar 




5(7+v^) 
1 Ar 


5(7-V7) 


2 


3 




1 


36 

A 
72 
A 


A^ 
144 
A2 

144 


-r Ar 

-■- 72 

1 Ar 

36 


3 







1 

2 


3.74329 

13.2216 

A 
31.0351 

A 


64 
A^ 

64 


v-^ 13.2216^"^-^ 
V-^ 3.74329 ''v 
V 13.2216^"^ 


Ar N 


31.0351'' 
Ar N 


31.0351 '' 
Ar \ 


3.74329'' 



5. Energy eigenvalues of a Hydrogenic Donor 

The Hainiltonian of an on-center impurity in a spherical quantum dot can be written 
in the effective-mass approximation as 



H 



^ w2 e2 
2m* tr 



(17) 



where m* is the effective mass and e is the dielectric constant of the material of the 
quantum dot. The donor is assumed to be at the center of the quantum dot of radius R 
with an infinite barrier height. This means that the wave function vanishes at r = i?. 
In atomic units, the radial Hamiltonian equation for the Coulomb Potential is given by 

2" 



rf2 /(/ + !) 



dr 



^(r) = E^(r), V(0) = V"!^) = 0. 



In this section, we use AIM in order to obtain the energy states of the confined atom. 
We note first that, Eq.Q now reads 

1 / + !' 



r(r) = 2 a + 



+ 



R-r 

{2l + 2)a-2 



fir) 



+ 



2/ + 2 
r{R — r 



2a 
R-r 



fir) 



(19) 
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Now, for a given R value, we have calculated the energy values a using ([! 
recursive relations (IHl) initiated with 



and the 



(20) 



Since the system is not exactly solvable, and since (5„(r) = Sn{a; r) = 0, we must choose 
a suitable value for ro < -R in order to initiate AIM. For our numerical results, we have 
fixed r = ro as ^ for R < 1 and 1 for R > 1. The corresponding results are shown 
in Tables 2 and 3. Comparison is made with results from ^7\, which are based on a 
variational method. 



' Ao = 2(a+^-^), 




„ _ /(2/+2)a-2 1 21+2 
\^ ^0 y r ' r(R-r) 


2a 
R- 



Table 2. Exact energy eigenvalues for the 1st state (in Rydbergs) found with AIM, 
for different values of R. Comparison is made with the results of Ref. [17]. 



R 


a 


Eaim 


EexactX^ 


0.1 


30.626 558 365 553 640 428 52i 


937.986 077 318 663 675 020 


937.986 


0.2 


14.904 352 306 411 647 880 03i 


222.139 717 673 638 207 696 


222.14(0) 


0.3 


9.653 226 131 767 727 449 50i 


93.184 774 751 043 322 516 


93.185(0) 


0.4 


7.019 089 585 470 665 068 m 


49.267 618 608 862 752 786 


49.268(0) 


0.5 


5.431 016 485 033 032 864 57i 


29.495 940 060 700 559 289 


29.496(0) 


0.6 


4.365 250 922 031 022 986 22i 


19.055 415 612 292 696 322 


19.055(0) 


0.7 


3.597 200 613 602 544 685 31i 


12.939 852 254 502 523 992 


12.940(0) 


0.8 


3.014 425 378 412 749 608 39i 


9.086 760 362 018 848 673 


9.0868(0) 


0.9 


2.554 286 411 599 184 224 80i 


6.524 379 072 480 237 168 


6.5244(0) 


1.0 


2.178 986 400 188 704 127 38i 


4.747 981 732 207 327 454 


4.7480(0) 


1.2 


1.593 307 889 331 964 142 70^ 


2.538 630 030 207 478 496 


2.5386(0) 


1.4 


1.137 633 609 176 779 056 96i 


1.294 210 228 728 584 474 


1.2942(0) 


1.6 


0.736 630 589 422 075 466 m 


0.542 624 625 272 314 320 


0.54262 


1.8 


0.255 171 521 889 990 724 15i 


0.065 112 505 583 654 015 


0.06511 


2.0 


0.5 


-0.25* 


-0.25 


2.2 


0.681 224 443 905 427 229 15 


-0.464 066 742 974 258 570 


-0.46407 


2.4 


0.782 812 885 873 172 513 99 


-0.612 796 014 289 084 615 


-0.61280 


2.6 


0.847 323 189 868 271 355 52 


-0.717 956 588 088 542 630 


-0.71796 


2.8 


0.890 692 535 878 441 133 67 


-0.793 333 193 469 568 146 


-0.79333 


3.0 


0.920 833 630 721 047 974 35 


-0.847 934 575 466 907 348 


-0.84793 


3.2 


0.942 234 555 027 191 412 81 


-0.887 805 956 687 289 402 


-0.88781 


3.4 


0.957 650 893 749 316 017 63 


-0.917 095 234 298 863 756 


-0.91710 


3.6 


0.968 866 519 413 519 397 41 


-0.938 702 332 440 467 559 


-0.93870 


3.8 


0.977 080 758 210 929 418 85 


-0.954 686 808 066 044 717 


-0.95469 


4.0 


0.983 122 883 548 157 499 83 


-0.966 530 604 156 044 052 


-0.96653 



Exact, see Table 1 
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Table 3. Exact eigenenergies (in Rydbergs) obtained with AIM for different values of 
R for the 2p state {£ — 1). Comparison is made with results of Ref. [IT], see also Ref. 



R 


a 


Eaim 


E exact [17] 


0.4 


10.811 856 798 907 243 508 22i 


116.896 247 440 076 786 588 403 


116.896 


1 


4.055 400 921 280 376 683 83i 


16.446 276 632 321 727 964 741 


16.446 


2 


1.775 397 862 793 768 780 65^ 


3.152 037 571 212 681 836 807 


3.1520 


4 


0.535 774 362 421 267 801 33i 


0.287 054 167 427 916 019 155 


0.28705 


8 


0.457 055 940 572 194 782 44 


-0.208 900 132 812 333 648 630 


-0.2089 



6. Critical cage radii 

In this section we report the results of our calculations for the critical cage radius, at 
which the total energy is zero, i.e., when the kinetic and potential energy contributions 
cancel one another. Such cage radii were first noted by Sommerfeld and Welker [5] in 
their detailed study on the variation of the binding energy of the Is state of the hydrogen 
atom, as a function of the sphere radius, R. Assuming that the surface of the spherical 
box is impenetrable, they showed that, as R decreases, the binding energy diminishes, 
and there is a critical value of the sphere radius at which the binding energy becomes 
zero. More systematic studies were later carried by Varshni [30] where the critical cage 
radius was first recorded; see also [31]. In this section we shall use AIM to calculate the 
cage radius at which the eigenvalues are zero. With ii^ = 0, we find, by using (jll), that 
a = 0; and then Eq. (TT9|) reduces to 

/ 1 ^ + l\ .,, . / 2 2£ + 2 



f"{r) 



fir) + 



fir) 



(21) 



R — r r j ' ^ ' ' \ r r{R — 

We now apply AIM with Aq and Sq extracted from (|2T]) . The termination condition 
6n{r) = 0, is then dependent (see (jH])) on both the variable r and the parameter R = re- 
Starting with r ^ -R/2, the iteration process converges quickly to the cage radius R ~ Vc- 
For example, if £ = and n = 5, AIM takes 26 iterations to yield Vc = 48.09774 initiated 
with starting value of r = 24. The calculated values of r^ found by using AIM are shown 
in Table 4 for different values of I and n. Similar tables can be easily constructed for 
arbitrary values of n and i. The numerical computation of eigenvalues and the critical 
cage radii in tables 2-4 were performed by using Maple version 10 running on an IBM 
architecture personal computer (Dell Dimension 4400). In many cases, we have removed 
some of the apparent divergence [37] experienced by AIM by increasing the number of 
significant digits that Maple uses in numerical computation; for the present results 
we have used Digits=50. In order to accelerate the computation we have written our 
code for the root-finding algorithm, instead of using the default procedure of Maple. 
All the numerical results reported using AIM in tables 2-4 are exact in sense that the 
numerical integration of the corresponding Schrondinger (confined) equation yields the 
same values. 
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Table 4. Exact values of the critical cage radius calculated by the use of AIM 
rc{AIM). Radii are in units of the Bohr radius. Comparison with |30| are also reported 



i 


n 


r^iAIM) 


re 


i 


n 


re{AIM) 


Tc 





1 


1.835 246 330 265 5 


1.8352 


1 


1 


5.088 308 227 275 


5.0883 




2 


6.152 307 040 211 8 


6.1523 




2 


11.909 696 568 004 6 


11.910 




3 


12.937 431 736 892 1 


12.937 




3 


21.174 431 228 262 4 


21.174 




4 


22.190 095 851 725 6 


22.190 




4 


32.900 106 781 876 


32.900 




5 


33.910 206 784 109 2 


33.910 




5 


47.090 674 929 020 9 


47.091 




6 


48.097 738 137 838 7 


48.098 




6 


63.747 459 484 409 4 


63.747 



7. Conclusion 



In this work, we have applied AIM to obtain the energy eigenstates of a confined 
hydrogen atom. The problem is similar to the confinement of electrons in a quantum 
dot. For certain cases, the eigenvalues under confinement are given by the roots 
of the Kummer (confiuent) hypergeometric functions whose analytic and numerical 
properties are well known |36] . The method can be easily adapted to the study of more 
highly excited states, without necessitating extensive further algebraic manipulation or 
numerical work [33] . We have also calculated the critical cage radii Tc at which various 
states become unbound. The numerical results are compared with the most accurate 
results in the literature. It is worth pointing that AIM provides a simple technique to 
obtain very accurate eigenenergies for a confined hydrogen atom, as well as the critical- 
cage radii, to any desired degree of precision. The method is easily realized by the use of 
any contemporary mathematical software. This makes the study of confined potentials 
more accessible. 
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